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Abstract 

We present a new analysis of the ratio e'/e within the Standard Model (SM) using a 
formalism that is manifestly independent of the values of leading (V — A) < 8 > {V — A) 
QCD penguin, and EW penguin hadronic matrix elements of the operators Q4, 
Q9, and Q 10 , and applies to the SM as well as extensions with the same operator 
structure. It is valid under the assumption that the SM exactly describes the data 
on CP-conserving K —>• 7T7 t amplitudes. As a result of this and the high precision 

now available for CKM and quark mass parameters, to high accuracy e'/e depends 

(1/2) (3/2) 

only on two non-perturbative parameters, B g 7 7 and B% 7 , and perturbatively 
calculable Wilson coefficients. Within the SM, we are separately able to determine 
the hadronic matrix element {Qa)q from CP-conserving data, significantly more pre¬ 
cisely than presently possible with lattice QCD. Employing B g 1 ^ = 0.57±0.19 and 
S< 3/2) = 0.76 ± 0.05, extracted from recent results by the RBC-UKQCD collabora¬ 
tion, we obtain e' je = (1.9 ±4.5) x 10~ 4 , substantially more precise than the recent 
RBC-UKQCD prediction and 2.9 a below the experimental value (16.6±2.3) x 10~ 4 , 
with the error being fully dominated by that on B g 1 ^. Even discarding lattice input 
completely, but employing the recently obtained bound Bq ~ < B$ ~ < 1 from 

the large-IV approach, the SM value is found more than 2 a below the experimental 
(1/2) (3/2) 

value. At Bq ~ = B$ ~ = 1, varying all other parameters within one sigma, we 
find e' /e = (8.6 ± 3.2) x 10 -4 . We present a detailed anatomy of the various SM 
uncertainties, including all sub-leading hadronic matrix elements, briefly comment¬ 
ing on the possibility of underestimated SM contributions as well as on the impact 
of our results on new physics models. 
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1 Introduction 


One of the important actors of the 1990s in particle physics was the ratio e'/e that 
measures the size of the direct CP violation in K L —>■ m r relative to the indirect CP 
violation described by £k■ In the Standard Model (SM), e' is governed by QCD penguins 
but receives also an important destructively interfering contribution from electroweak 
penguins that is generally much more sensitive to new physics (NP) than the QCD penguin 
contribution. Reviews on s' je can be found in 0R 

A long-standing challenge in making predictions for e'/e within the SM and its ex¬ 
tensions has been the strong interplay of QCD penguin contributions and electroweak 
penguin contributions to this ratio. In the SM, QCD penguins give a positive contribu¬ 
tion and electroweak penguins a negative one. In order to obtain a useful prediction for 
e'/e, the relevant contributions of the QCD penguin and electroweak penguin operators 
must be know accurately. 

As far as short-distance contributions (Wilson coefficients of QCD and electroweak 
penguin operators) are concerned, they have been known already for more than twenty 
years at the NLO level [6-11 and present technology could extend them to the NNLO 


level if necessary. First steps in this direction have been taken in 12-14 


The situation with hadronic matrix elements is another story and even if significant 
progress on their evaluation has been made over the last 25 years, the present status 
is clearly not satisfactory as we will discuss below. But, already in 1993, an approach 


has been proposed in 10 which, as far as e'/e is concerned, avoids direct calculation of 
some of the most difficult hadronic matrix elements. It assumes that the real parts of 
the isospin amplitudes A 0 and A 2 , which exhibit the A/ = 1/2 rule, are fully described 
by SM dynamics and their experimental values are used to determine to a very good 
approximation hadronic matrix elements of all (V — A) <g) {V — A) operators, among them 
the so-called QCD penguin operator. While not as important as the (V — A)<S> (V + A) 
QCD penguin and electroweak penguin operators, Qe and Q §, the operator Q± has been 
known since the early days of analyses of e'/e |9,10 15 to be responsible for a significant 
part of the suppression of this ratio. In the presence of a partial cancellation of the 
positive contribution of Q§ to e'/e by the one of Q$, an accurate determination of the 
contribution from Q 4 and from the electroweak penguin operators Qg and Q io to e'/e by 
means of CP-conserving data was an important virtue of our approach. 

Another virtue of our approach is based on the fact that in the SM the amplitudes 
ReA 0 an d ReA 2 originate already at tree-level. Similar to the observables used for tree- 
level determination of CKM parameters, also relevant for e'/e, they are expected to be 
only marginally affected by NP contributions. Whether NP could contribute to ReAo and 
ReA 2 at some level is an interesting question, to which we will return briefly in Section 6. 
But, for the time being we assume that they are fully dominated by SM dynamics. 

With the contribution of (V — A) <g) (V — A) operators being determined from the 
data on ReA 0 and ReA 2 it was possible to write down an analytic formula for e'/e that 
incorporated all NLO QCD and QED corrections and summarised the remaining dominant 
hadronic uncertainty in terms of two parameters B g 1 ^ 2 '* and B g 3 ^ that parametrise the 
relevant matrix elements of the dominant operators Qq and Qs and have to be calculated 


using a non-perturbative framework like lattice QCD or the large-A approach 16,17 . 
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They cannot be extracted from CP-conserving data as their contributions to ReAo and 
RehL 2 are marginal at n ~ m c used in the approach of 10 . In fact one of the reasons for 


choosing the value /i = m c was to eliminate them from the determination of the matrix 
elements of (V — A) ® (V — A) operators from the CP-conserving data. 



One new aspect of the present paper is the realisation that under the assumption 
that NP contributions to ReAo and ReA 2 are negligible, the leading contributions of 
(V — A) 0 (V — A) operators to s' je can be entirely expressed in terms of their Wilson 
coefficients. Furthermore, we derive a formula for s' je which under the above assumption 
can be used in any extension of the SM in which the operator structure is the same as 
in the SM. NP enters only through the modified values of the Wilson coefficients and the 
dominant non-perturbative uncertainties are contained in 


B, 


( 1 / 2 ) 


B 


( 3 / 2 ) _ z +(f i )(Q+(f 1 )) 0 

’ q - „ ' 


( 1 ) 


The ratio q, involving matrix elements of current-current operators Q± and their Wilson 
coefficients z±, enters the determination of the contribution of (V — A)<S>(V — A) operators 
from CP-conserving data and its range will be estimated in Section 2. But for 0 < q < 0.1 
obtained from QCD lattice and large- N approaches the dependence of s' je on q is very 
weak. 

As far as the parameters B g 1 ^ and B g 3 ^ are concerned, B g 1 ^ 2 ' = B^ 2) = 1 in 
the large- N limit of QCD. The study of 1/N corrections to the large- N limit indicated 
that B § 7 ; is suppressed below unity [22|, but no clear-cut conclusion has been reached 
in that paper on B^ 2 \ Moreover, the precise amount of suppression of B g 3 ^ could 
not be calculated in this approach. Fortunately, in the meantime significant progress 
has been achieved in the case of the matrix element (Q %)2 by the RBC-UKQCD lattice 
, which allowed to determine B g 3 ^ 2 '* to be 


collaboration 


23 


21 


B 


( 3 / 2 ) 


(m c ) = 0.76 ± 0.05 


in agreement with large -N expectations 22,24 


(RBC-UKQCD), 
but with higher precision. 


( 2 ) 


But also some progress on B g 1 ^ 2 "* has been made, both by lattice QCD and the large- 
N approach. In particular, very recently the RBC-UKQCD lattice collaboration |25 


presented their first result for the matrix element (Qq )o from which one can extract (see 
below and |24] ) 

Rg 1/2) (m c ) = 0.57 ± 0.19 (RBC-UKQCD). (3) 

This low value of B^ 2 ' 1 is at first sight surprising and as it is based on a numerical 
simulation one could wonder whether it is the result of a statistical fluctuation. But the 
very recent analysis in the large-lV approach in 24 gives strong support to the values 
in (2) and (3). In fact, in this analytic approach one can demonstrate explicitly the 
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suppression 


of both B^ 2 ' 1 and B^ 2 ' 1 below their large-iV limit B^ 2) = B^ 2) = 1 and 


derive a conservative upper bound on both B g 1 ^ and B^ 2 ^ which reads 

(large- N). 


24 


B^ /2) < Bf /2) < 1 


(4) 


While one finds B^ 2 \m c ) = 0.80 ±0.10, the result for B^ 2 ^ is less precise but there is a 
strong indication that B g 1 ^ < B^ 2 ' 1 in agreement with (2) and (3). For further details, 
24 and Section 3; below. 


see 


Employing the lattice results of (2) and (3]), in our numerical analysis we find 

e'/e = (1.9 ±4.5) x 10~ 4 , (5) 

consistent with, but significantly more precise than the result obtained recently by the 
RBC-UKQCD lattice collaboration [25], 


(e'/e) S M = (1-4 ± 7.0) x 10 


-4 


( 6 ) 


This is even more noteworthy considering the fact that our result comprises also uncertain¬ 
ties from isospin corrections and CKM parameters which were not considered in the error 
estimate of 25 . Our result differs with close to 3 cr significance from the experimental 

collaborations, 

(7) 


world average from NA48 |26| and KTeV |27[ 

(e7 £ )exp = (16.6 ± 2.3) x 10" 4 


suggesting evidence for new physics in K decays. 

But even discarding the lattice results, varying all input parameters, we find at the 
bound B { 6 1/2) = B£ /2) = 1, 


(e'/e)sM = (8.6 ± 3.2) x 1(T 4 , 


( 8 ) 


still 2 cr below the experimental data. We consider this bound conservative since employing 
the lattice value in (2) and B^ 2 ^ = B^ 2) = 0.76, instead of (8), one obtains (6.0 ±2.4) x 
10 -4 . 

(' 3 / 2 ) 

This already shows that with the rather precise value of B$ 7 ; from lattice QCD, the 
final result for s'/e dominantly depends on the value of B g 1 ^ and both lattice QCD 

indicate that the SM value of s' je is significantly be 


24 


and the large-iV approach 
the data. 

The two main goals of the present paper are: 


25 


ow 


• Derivation of a new version of our formula for e'/e which could also be used beyond 
the SM and which appears to be more useful than its variants presented by us in 
the past. 

• Demonstration that our approach provides a substantially more accurate prediction 
for e'/e in the SM than it is presently possible within lattice QCD and that the 
upper bound in (8) is rather conservative. 
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It should be stressed that assuming dominance of SM dynamics in CP-conserving data, 
our determination of the contributions of (V — A) <g) (V — A) operators to e'/e is basically 
independent of the non-perturbative approach used. The RBC-UKQCD lattice collab¬ 
oration calculates these contributions directly and we will indeed identify a significant 
difference between their estimate of the Q 4 contribution to e'/e and ours. 

Our paper is organised as follows. In Section 2, we derive the analytic formula for e'/e 
in question using the strategy of 10 but improving on it. Using this formula, we present 
a new analysis of e'/e within the SM exhibiting its sensitivity to the precise value of B2^ 
and the weak dependence on q. In Section |3j we perform the anatomy of uncertainties 


iction of e'/e in the SM, including a discussion of 
we extract from the lattice-QCD results of 25 the 


affecting e'/e and present the prec 
its B g 1 ^ dependence. In Section 
values of the most important hadronic matrix elements and compare them with ours. This 
allows us to identify the main origin of the difference between (J5| and In particular, 
we point out an approximate correlation between the contribution of the Q 4 operator 
to e'/e and the value of ReAo valid in any non-perturbative approach. In Section [5j, 
we investigate if thus far neglected SM contributions could bring our result for e'/e into 
agreement with the experimental findings. A brief general discussion of the impact of 
possible NP contributions to ReAo ,2 and ImA 0 ,2 and of the implications of our results for 
NP models is given in Section 6. The summary of our observations and an outlook are 
presented in Section [7| I 11 Appendix [A] we discuss the sub-leading contributions to our 
prediction for e'/e and in Appendix B, for completeness, an updated analytic formula for 
e'/e in the SM is presented in the form used in several of our papers in the past (e.g. j|2l] ) 
that is equivalent to the one derived in Section 2, but exhibits the m t , a s , m s and rn c i 
dependences more explicitly. 


2 Basic formulae 

2.1 Effective Hamiltonian 


H eS = C ^V ud V* s ^(z i (n)+Ty i (n))Q i (fi), t = — vywry ■ (9) 

The contributing operators are given as follows: 

Current—Current: 



Q 1 — (s a Up)v-A (Upd a )v- 

-A Q 2 ~ {su)v-A ( ud) v _ A 

(10) 

QCD Penguins: 



Q 3 — 

( sd) V -A ^2 (QQ)v-a 
q=u,d,s,c,b 

Q 4 = ( S a dp)v-A '22 ( QpQa)v-A 

q=u,d,s,c,b 

(11) 

Q 5 = 

(sd) V - A (m)v+a 

q=u,d,s,c,b 

Qe = (, s a dp)v-A 22 (QpQo)v+a 

q=u,d,s,c,b 

(12) 


We use the effective Hamiltonian for AS = 1 transitions of [6Ml 
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Electroweak Penguins: 


Qy = a (■ sd)v-A e i ( 99 ) 


V+A 


q=u,d,s,c,b 


Q 9 = ^ {sd)v-A ^2 e i (vtiv-A 

q=u,d,s,c,b 


Qs 2 ^adp)v—A / v 

q=u,d,s,c,b 


Q 10 — A ( S Q dp)v-A ^2 


e q ( qpq a )v+A (13) 


a V-A 


(14) 


q=u,d,s,c,b 

Here, a, (3 denote colour indices and e q denotes the electric quark charges reflecting the 
electroweak origin of Q 7 ,..., Q w . Finally, ( sd)y-A = s a 7 M (l — 7 s)d Q . 

The Wilson coefficients z t and i/i have been calculated at the NLO level more than 


twenty years ago 10,11 , and some pieces of NNLO corrections are also available 12114 


In Table [TJ we collect values for z\^. and yt at n — m c , used in our approach, for three 
values of oi s (M z ) and m t = 163 GeV, in the NDR-MS scheme. 



a s {M z ) =0.1179 

a s (Mz) = 0.1185 

a s (M z ) = 0.1191 

Z\ 

-0.4036 

-0.4092 

-0.4150 

Z2 

1.2084 

1.2120 

1.2157 

2/3 

0.0275 

0.0280 

0.0285 

2/4 

-0.0555 

-0.0563 

-0.0571 

2/5 

0.0054 

0.0052 

0.0050 

2/6 

-0.0849 

-0.0867 

-0.0887 

2 / 7 /« 

-0.0404 

-0.0403 

-0.0402 

2/s/« 

0.1207 

0.1234 

0.1261 

2 / 9 / a 

-1.3936 

-1.3981 

-1.4027 

2/io/a 

0.4997 

0.5071 

0.5146 


Table 1 : AS = 1 Wilson coefficients at fx = m c — 1.3 GeV for three values of a s (M z ) and 
m t = 163 GeV in the NDR-MS scheme. 


2.2 Basic formula for s'/s 

Our starting expression is formula (8.16) of 29 which we recall here in our notation 1 


£ 


U+ 


y/2 


£k\ . 


Imy4 0 ft n \ Im ^2 
ft — “effj — 


ReHn 


Re/lo 


where 129 


ijj + = a 


Re/D 

Re/ln 


= (4.53 ±0.02) x 10' 2 , a = 1.017, fi eff = (6.0 ± 7.7) x 10 


-2 


(15) 


(16) 


Here a and O e ff summarise isospin breaking corrections and include strong isospin violation 
(rn u 7 ^ rrirf), the correction to the isospin limit coming from A I = 5/2 transitions and 

Hn order to simplify the notation we denote Refe'/e) simply by e'/e, which is real to an excellent 
approximation. The latter is a model-independent consequence of the experimentally known values of 
the (strong) phases of e’ and e. 
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electromagnetic corrections 29,30 . The amplitudes ReAo ,2 are then extracted from the 


branching ratios on K —> tttt decays in the isospin limit. Their values are given in (47) 
below. In the limit a = 1 and h2 e ff = 0 formula (15) reduces to the one used in 25 , where 


all isospin breaking corrections except electroweak penguin contributions have been set 
to zero. 

The quantity fi e ff includes, in addition to other isospin breaking corrections, elec¬ 
troweak penguin contributions that are then not included in ImAo- Here we prefer to 
include these contributions to ImAo and therefore, instructed by the authors of 29 , we 
remove them from n e g. However, we keep in n e g- their term A 0 in the limit of a = 0. 


Using Table 4 of 29 , we then obtain the modified 


n eS = (14.8 ± 8.0) x 10 


-2 


(17) 


As the second term in (15) is an isospin breaking effect by itself, strictly speaking in this 


term the parameter a should be set to unity if we want to remove higher order isospin 
braking corrections. In addition, in order to remove the effects of fi e g- ^ 0 in electroweak 
penguin contributions to ImA 0 , we write 


ImA 0 = (Inn4 0 ) QCDP + b (ImA 0 ) 


EWP 


b = 


o, (i ii e g) 


(18) 


with the first term including the contributions from Q 3-6 and the second from Q 7 _io- 
Except for the tiny corrections due to a ^ 1, this procedure is equivalent to multiplying 
the coefficients t/ 3 _ 6 by (1 — H e ff) leaving r/ 7 _ 10 unchanged. 

Our final basic formula which we will use in what follows then reads 


UJ+ 


V2\e K \ 


ImA 0 

ReAn 


(1 — H e ff)- 


1 ImAo 


a ReA 2 


(19) 


with (a; + ,a), fl e ff and ImA 0 given in (16), (17) and (18), respectively. lmA 2 contains only 
contributions of the electroweak penguin operators Qt-io- 

The crucial theory task for a precision SM prediction is to determine the real and 
imaginary parts of the (strong-)isospin amplitudes 


A/ = {{'Kn)i\'H e fi\K) 


( 20 ) 


entering (19) in terms of the Wilson coefficients and hadronic matrix elements of the 


operators in the weak Hamiltonian (9). 

2.3 Hadronic matrix elements 


The hadronic matrix elements of the operators Qi entering the isospin amplitudes, 


{Qi)i = ((tttt)/ \Qi\ K ), 7 = 0,2, (21) 

generally depend on the scale fi and on the renormalisation scheme used for the operators. 
These two dependencies are cancelled by those present in the coefficients Ci(fi) so that the 
effective Hamiltonian and the resulting amplitudes do not depend on /i and on the scheme 
used to renormalise the operators. We will work exclusively in the NDR-MS scheme and 
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10 also extensive discussion of scales above m r can be 


for scales /i < m c , although in 
found. 

For /j, < m c , when the charm quark has been integrated out, only seven of the operators 
listed above are independent of each other. Eliminating then Q 4l Qg and Q 10 in terms of 
the remaining seven operators results in the following important relations in the isospin 
limit 


10 


(^4)0 — (^3)0 + {Q 2)0 — (Qi)o, ( 22 ) 

(^ 9)0 = 2 (Qi)o — 2 ^ 3 )° ’ (23) 

(Qio)o = (<^ 2)0 + 2 (Qi)o — 2 (^ 3 ) 0 ’ (24) 

(^9)2 = (Q 10)2 = 2 ($1)2 j ( 25 ) 

where we have employed 

(Q 1)2 = (<32)2 • ( 26 ) 


As stressed in 10 , in the NDR-MS scheme the relation (22) receives an 0(a s ) cor¬ 


rection due to the presence of evanescent operators which have to be taken into account 
when using Fierz identities in its derivation. The other relations above do not receive 
such corrections. The complete expression for (< 54)0 in the NDR-MS scheme reads 


10 


(Q4)o — (^3)0 + (^2)0 — (Qi)o — 4" ($4)0 “ g (Q3)o ~ g (Qs)o) , (27) 


which of course then has to be solved for (Q 4 ) 0 . However, due to the partial cancellation 
between the matrix elements (Q 4 )0 and (Qe) 0 , and the smallness of the matrix elements 
of Q 3 and Q 5 , this correction affects the determination of (Q 4 ) 0 by at most few percent 
and can be neglected. This procedure is supported both by the results on hadronic matrix 

Setting the contribution of Q 3 to zerc 2 and using the operators 


elements RBC-UKQCD collaboration 25 and the large -N approach 24 


the formulae (22)-(25) read 


Q ± = 2 {Qz 4 = Qi ), 

(28) 

(<34)0 = 2 ( Q ~) 0, 

(29) 

(^9)0 = 2 {( Q +) 0 — ( Q -) 0) j 

(30) 

3 1 

(Qio)o = 2 (Q+)o + 2 ( Q-)o , 

(31) 

(<^9)2 = (Q 10)2 = 2 (Q+) 2 > 

(32) 


2 In our numerical analysis below, all operators will be taken into account. 
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which reduces the number of independent (V — A) (g) (V — A) matrix elements entering 
Rehlo ,2 and ImA^^ to three. On the other hand, to an excellent approximation the am¬ 
plitudes RehL 0 and ReA 2 at fi = m c are fully described by the operators Q- and Q + , so 
that we can write 


ReA 0 — — y= V ud V* s (z+{Q +) 0 + z_(QAo) , 

(33) 

ReA 2 = C ^V ud V: s z + (Q + ) 2 - 

(34) 

Introducing the ratio 


Q — ( \ //n ( \\ j *2 i *1 i 

z-W(Q-W) 0 

(35) 


allows us to express the ratios involving only (V — A) ® (V — A) operators that will enter 
our basic formula for e'/e as follows: 


f Im ^o \ = Imr [% 4 ~ b(3yg - y w )] fj 3g(y 9 + 2/io) 
yRe^o/v-A 2(1 + q)z^ 2(1 + q)z + 


(36) 


f Im A 2 \ _ 3(y 9 + y io) 

\ReA 2 ) v _ A 2z+ 


(37) 


Besides the CKM ratio r, the first ratio depends only on Wilson coefficients and the single 
hadronic ratio q to which we will return below. On the other hand the second ratio is free 
from hadronic uncertainties, being fully determined by the Wilson coefficients z + , yg, yio 
and by r. 

The remaining contributions to Inn4 0 and ImA 2 are due to (V — A) 0 (V + A) operators 
and are dominated by the operators Q 6 and Q 8 , respectively. We find this time 


f Im/If 


\ReA 0 
/ Im A 2 


Gf t n (Qg )o 

= - ~7= ImA « 2/6 

6 v 2 ReA 0 


Gf t , e 
= ~^ lmXty s 


.eff 


(Qs)* 


(38) 


(39) 


\ReA 2 J 8 ^2 ™ ReA 2 

Contributions from Q 3 and Q 5 are very suppressed but can and have been included in 
our numerical error estimate. (See Appendix 0 ) We have also taken into account the 
small effect of (< 37 ) 2 , for which a relatively precise lattice prediction exists 23 , through 
the substitution 

2/8 -t 


pff 

2/8 =V8+P72y7 


(40) 

Here = ((^ 7 ) 2 /(QA 2 = 0.222 for central values 
(In our numerics, we have added the corresponding errors linearly and attribute 
a 15% uncertainty to this contribution.) 

The matrix elements of the Qq and Q 8 operators are conveniently parametrised by 


which is included in writing (39) 
of [23 


(Qe (/■*)> 0 = -4 h 


m K 


m s (p) + m d (p) 


(Fi 


I< 


FAB , 


( 1 / 2 ) 


(41) 


(Q8(h))2 


V2h 


m 2 K 


m s (p) + m d (n)_ 


Fir B 


( 3 / 2 ) 


(42) 
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with 31,32 


B, 


( 1 / 2 ) 


= B, 


( 3 / 2 ) 


= 1 


(43) 

in the large-A limit. As had been demonstrated in 10 , B g 1 ^ and B^ 2 ^ exhibit a 
very weak scale dependence. The dimensionful parameters entering (41), (42) are given 


by 33,34 


m K = 497.614 MeV, 


= 130.41(20) MeV, ^ = 1.194(5), 

Ftt 


m s (m c ) = 109.1(2.8) MeV, m d {m c ) = 5.44(19) MeV , 


(44) 

(45) 


In 34 , the light quark masses are presented at a scale of 2 GeV, and we have evolved 


them to /i = m c = 1.3 GeV with the help of the renormalisation group equation. For the 
comparison with lattice results below, we also need their values at /i = 1.53 GeV, which 
are found to be 


m s (1.53 GeV) = 102.3(2.7) MeV, m d (1.53 GeV) = 5.10(17) MeV , 
Below, we will neglect the tiny errors on mjf, F^, and F n . 


(46) 



ReA 0 = 33.22(1) x l(T 8 GeV, 
which display the A/ = 1/2 rule 


ReA 2 = 1.479(3) x 1(T 8 GeV, 


(47) 


ReA 0 1 
—A = - = 22.46. 
ReA 2 w 


(48) 


10 


We also note that while equation (41) is identical to (5.10) in 10 , the definition of B 
in the present paper differs from 


>( 3 / 2 ) 


[cf (5.18) there]. This is to ensure that B^ 2) = 1 

/o/o) 

and B% =1 both correctly reproduce the large-A limit of QCD. In contrast, (5.18) 


m 


10 was based on the so-called vacuum insertion approximation, in which additional 

Such terms misrepresent the large-A limit 


( 3 / 2 ) 


terms appear in the normalisation of B 
of QCD. With our conventions, 1/A corrections in (41) and (42) are represented by the 
departure of Bf } l,/2> and B^ 2 ^ from unity. They have been investigated in 
24 with the result summarised in (4). 


recently in 
details. 


22 


and very 


We refer to this paper for further 


We now turn to the parameter q which enters (36). We first note that, like B q 1 ^ and 
B g ' , it is nearly renormalisation-scale independent. Its value can be estimated in the 
large-A approach 17 ; as this approach correctly accounts for the bulk of the experimental 
value of ReAo ; the ensuing estimate can be considered a plausible one. In the large-A 
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limit, corresponding to ji = 0, one finds first (Q + (0))o/(Q-(0))o = 1/3. Using the meson 
evolution in (TF) up to p = 1.0 GeV and then quark evolution up to p = m c , multiplying 
the result by z + (m c )/ Z-(m c ), we obtain q ~ 0.1. On the other hand the results of the 
RBC-UKQCD collaboration [25] are consistent with a value of zero ( q = 0.029 ± 0.087). 
As the large -N approach gives ReA 0 below the data while 25 above it, we expect the 


true value of q at p = m c to lie between these two estimates and will take q in the range 

0<g<0.1. (49) 


We consider this a credible range, but already mention that our phenomenological results 
below would change very little even if we enlarged this range by a factor of a few: q is 
simply too small to introduce a large error on e' je. 

Our input parameters including sub-leading hadronic parameters defined in Appendix 
[Aj are collected in Table |2. Regarding IrnA t , we choose a central value between the UT- 
fit [3§ and CKMfitter [37] determinations and an error slightly larger than that obtained 
from either fit. This is to account for the very small errors on V ud and V us , which we fix to 
PDG central values 33]. The Wilson coefficients in Table |l] come with an additional un¬ 


certainty from unknown higher-order corrections. In particular the threshold corrections 
at m c can be substantial even at NNLO. This can for example be seen in the perturbative 
convergence of £x 14,38]. We use a scale variation to establish the typical size of higher 
order corrections and estimate a 10% uncertainty for each Wilson coefficient yz - pio of 
Tabled 



value range 

comment 

R ( l/ 2 ) 
n 6 

R ( 3/2) 
n 8 

q 

R (l/ 2 ) 

n 8 

P72 

P3 

P5 

P70 

0.57 ±0.19 

0.76 ±0.05 

0.05 ±0.05 

1.0 ±0.2 

0.222 ±0.033 

0 ± 0.5 

0 ± 0.5 

0 ± 1/3 

Eq. (1 
Eq. | 

Eq. (4 

1 ) and surroi 
1 ) and surroi 
see (35 
defined in 
0 ) and surro 
see Appc 
see Appc 
see Appc 

111 c 

111 c 

M 

Ec 

un< 

me 

mc 

mc 

in; 

in; 

49 

f ( 
iin 

ix 

ix 

ix 

?c 

? c 

87 

g 

A 

A 

A 

iiscussion 

iiscussion 

) 

discussion 

ImAi 

m t (m t ) 

m s (m c ) 

m d (m c ) 

s 2 

b W 

(1.4 ±0.1) x 10~ 4 

(163 ± 3) GeV 
(109.1 ±2.8) GeV 
(5.4 ± 1.9) GeV 
0.1185 ±0.0006 

0.23126 

see t 

calculated from po 
value from 
value from 
from 

MS scheme va 

ex 

1 c 

33 

33 

3e 

luc 

nass value 33 
evolved 
evolved 

from 33 

^efF 

(14.8 ±8.0) x 10" 2 

from 

21 

) 

V3 - VlO 

Vi x (1±0.1) 

see Text 


Table 2: Input parameter ranges, grouped into: hadronic matrix elements, parametric, 
isospin breaking and NNLO. The (numerically unimportant) ratios £>72, Pz, P5, P70 are 
defined in Appendix [AJ . The remaining parameters (iv, TV, mjf, V u a, V us , « em , Gf , £k) 
are fixed at their central values. 
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2.4 Convenient formula for s'/s 


Before turning to quantitative phenomenology, in order to make easier connection 
with the phenomenological literature and aid discussion of our results, we summarise 
the discussion so far in a concise formula (derived first in [To] ) for e'/e that exhibits 
the sensitivity to the two most important hadronic matrix elements B g / 7 and B g 7 7 
transparently. 

Using the effective Hamiltonian (9} and the experimental data for ui, Re+o and £k, 
we hnd 


£ 


lmA t • 


a(l — O eff ) p(l/ 2 )_p( 3 / 2 ) 


(50) 


where 


pd/2) = Y J P i' m = rJ2m(Qi)o, (51) 

P m) = Y. = - E . (52) 

with 

r = — -——— . 53 

2 \ek \ ReH 0 

In ( fHTj ) and ([52) the sums run over all contributing operators. Therefore in pd/ 2 ) j n the 
case of EWP contributions we have to take into account the correction 6^1 defined in 

Writing then 


pw = 4 1/2) + 4 1/2) p^ 1/2) , 


(54) 


p(3/2) _ J3/2) (3/2) r (3/2) 

r — u Q -f u g -Do , 


(55) 


with the parameters Pq 1//2 ^ and B^ 2 ' 1 taken at /i = m c and using the expressions (36)-(42) 
we hnd: 


( 1 / 2 ) 

a 0 = r l 


[4// 4 - 6(32/9 - pro)] b 3< 7 ( 2/9 + 2 / 10 ) 


2(1 + q) Z - 

{Qq) o 


2(1 + q)z+ _ 


+ r 2 by s 


{Qs)o 
ReHn ’ 


( 1 / 2 ) _ 

06 ~~ r2 2/6 Pg 1/2) ReH 0 

(3/2) _ 3(2/9 + 2/io) 

0 " fl 2z + ’ 

,.(3/2) _ _ „eff W 2 

8 “ 298 B ( 8 «Re4 2 


where 


cu 


cu 


ri = 


V^|£at| ’ 


r 2 


- 


■g f , 


(56) 

(57) 

(58) 

(59) 

(60) 
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and {Qq) o, (Qs)o, and (Q 3)2 are given in (41), (87) and (42), respectively. The second 
term in (56) proportional to q amounts at most to a 2% correction and could be safely 
neglected. is defined in (40). a and a ^ 2 ' 1 receive further small corrections which 
can be extracted from the expressions in Appendix [A] Apart from that, the coefficients 
a- 1 ^ and af^ depend only on q, a s , m t , and the renormalisation scheme considered. 
The dependencies on a s and m t are given in the NDR-MS scheme in Table p|> 


ot s (.M z ) 

m t [GeV] 

fl d /2 ) 

u o 


a (3/2) 

u 0 

»8 


160 

-2.93(12) 

17.23 

-0.82 

6.96 

0.1179 

163 

-2.90(12) 

17.25 

-0.84 

7.27 


166 

-2.87(12) 

17.26 

-0.85 

7.58 


160 

-2.95(12) 

17.61 

-0.82 

7.13 

0.1185 

163 

-2.92(12) 

17.63 

-0.84 

7.44 


166 

-2.89(12) 

17.64 

-0.85 

7.76 


160 

-2.98(12) 

18.00 

-0.82 

7.31 

0.1191 

163 

-2.95(12) 

18.02 

-0.84 

7.62 


166 

-2.92(12) 

18.03 

-0.85 

7.95 


Table 3: The coefficients a and af^ 2) in the NDR-MS scheme for different values of 
a s (Mz) and m t . The uncertainty shown for a ^ 2 ' 1 only includes the variation of q. 

I11 summary the ratio e'/e is governed by the following four contributions: 

i) The contribution of {V — A) ® (V — A) operators to pd/ 2 ) j s represented by the 


11 


first term in (54). As seen in (56) this term is governed by the operator Q 4 and 


includes also small contributions from {V — A) ® {V — A) electroweak penguin 
operators. We find that this term is negative and only weakly dependent on q. 
Also the dependences on a s and renormalisation scheme (see jTo] ) are weak. These 
weak dependences originate from the fact that in our approach the matrix elements 
entering the first term in Pd/ 2 ) cancel out. The weak dependence on m t results from 
the contributions of sub-leading electroweak penguin operators and is exhibited in 
the formulae in Appendix Bj As pointed out in 10 , the suppression of e'/e through 
4 1/2) increases with increasing ReA 0 , a feature which in the next section will help 
us to partly understand the result in (6). 

The contribution of (V — A) ® (V + A) QCD penguin operators to pd/ 2 ) j s given by 


the second term in (54). This contribution is large and positive and is dominated 


hi) 


by the operator Q e . The coefficient a ^ 2 ' 1 depends sensitively on a s , but as in the 
last two decades the precision on a s increased, this uncertainty is small in 2015 as 
can be seen from Table j3i 

The contribution of the (V — A)®(V — A) electroweak penguin operators Q 9 and Q\o 
to pd/ 2 ) is represented by the first term in pd/ 2 ). As in the case of the contribution 
i), the matrix elements contributing to a q 3 ^ cancel out in the SM. Consequently, 
the scheme and a s dependences of Oq ~ are weak. As seen in (58) the sizable rri t - 

( 3 / 2 ) 

dependence of ag ' results from the corresponding dependence of 1/9 + 1/10 but again 
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( 3 / 2 ) 

the precision on m t increased by much in the last two decades, a q 7 7 contributes 
positively to e'/e. 


iv 


The contribution of the (V — A)®(V + A) electroweak penguin operators Q 7 and Q 8 
to pAA) i s represented by the second term in (55). This contribution is dominated 
by Qs and depends sensitively on m t and a s . It contributes negatively to e'/e. 


The competition between these four contributions is the reason why it is difficult to 
predict e'/e precisely. In this context, one should appreciate the virtue of our approach: 
the contributions i) and iii) can be determined rather precisely by CP-conserving data so 
that the dominant uncertainty in our approach in predicting e'/e resides in the values of 


B, 


( 1 / 2 ) 


and B 


(3/2) 


3 Prediction for e'/e in the SM 


3.1 Prediction for e'/e and discussion 


We begin our analysis by employing the lattice values in (2) and (3). Varying all parame¬ 
ters within their input ranges and combining the resulting variations in e'/e in quadrature, 
we obtain: 


(^)sm = (1.9±4.5)x10- 4 . 

Comparing to the experimental result (£ / /e) eX p = (16.6±2.3) x 10 -4 (average of NA48 


and KTeV 27 28)), we observe a discrepancy of 2.9 a significance. 



quantity 

error on e'/e 

quantity 

error on e'/e 

d(1 /2) 

n 6 

4.1 

m d {m c ) 

0.2 

NNLO 

1.6 

q 

0.2 

^eff 

0.7 

I/ 2 ) 

0.1 

P3 

0.6 

IniA t 

0.1 

o(3/2) 

n 8 

0.5 

P72 

0.1 

P5 

0.4 

P70 

0.1 

m s (m c ) 

0.3 

ot s (M z ) 

0.1 

m t (m t ) 

0.3 




Table 4: Error budget, ordered from most important to least important. Each line shows 
the variation from the central value of our e'/e prediction, in units of 1CT 4 , as the cor¬ 
responding parameter is varied within its input range, all others held at central values. 

A detailed error budget is given in Table [4} It is evident that the error is dominated 
by the hadronic parameter Bq 2 \ Uncertainties from higher-order corrections are still 
significant yet small if compared to the deviation from the experimental value. All other 
individual errors are below 10” 4 , with the third most important uncertainty coming from 
the isospin breaking parameter at a level of 0.7 x 10~ 4 and about six times smaller 
than the error due to B/^^. If matrix elements are taken from a lattice calculation, the 
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m s dependence is only an artifact of our parametrisation in terms of B^ 2 ^ and B^ 2 \ 
Therefore including the m s variation in our error estimate for e'/e leads to a slight, but 
negligible, overestimate of the total error. At the same time, the small m s dependence 
we do find in the final result shows that this is no longer a relevant source of uncertainty 
in non-lattice approaches (like the large-iV approach in particular) in which B^ 2 ^ and 
Bg i, ' 2 ' > are directly calculated. 

At this stage it is important to emphasise that the results for B\ 7 7 and Bq 7 in 
(2) and |3) receive strong support from the large-A 7 approach as recently demonstrated 
in 24 . In particular the smallness of the matrix element (Qe )o with respect to (Qs )2 is 
the result of the chiral suppression of (Q e ) 0 , signalled by F K — F n in (41). As seen in (42) 


no such suppression is present in (<3 8 ) 2 . But in addition it is possible to demonstrate that 
both B^ 2 " 1 and B g 3 ^ are below unity as given in (4). Moreover, while B^ 2] = 0.8 ±0.1 
is found in this approach, the values of B g 1 ^ 2 '* are in the ballpark of the lattice result 
and consequently give a strong support for B$ 7 7 < B 8 7 " 7 as indicated by the lattice 
data. But as present calculations by lattice QCD and in 24 are not precise enough, at 
this moment, we cannot exclude that B^ 2 ' 1 could be as large as B 8 3 ^ and this leads 
conservatively to the bound in (4). 

For these reasons it is instructive to consider other values of the parameters B^ 2 ' 1 and 

(3/2) 

B g 7 7 than those obtained by RBC-UKQCD collaboration which are, however, consistent 
with the large- N bound in (4). Of particular interest is the choice B t 


( 1 / 2 ) 


= B 


(3/2) 


= 1 


which corresponds to the saturation of this bound and the choice in which the bound on 
(1/2) • ( 3 / 2 ) • • • n 

B g ~ is saturated when B § 7 7 is fixed to the central lattice value in (2). Using the same 

input for the remaining parameters, we find 


(eVe)sM = (8-6 ± 3.2) x HT 4 , (£< 1/2) 


= ^ 3/2) = 1 ), 


(eVe)sM = (6.0 ± 2.4) x 10 


-4 


( 5 , 


( 1 / 2 ) 


= B, 


(3/2) 


= 0.76). 


(62) 
(63) 

We observe that even for these values of B^ 2 ' 1 and B g 3 ^ the SM predictions for e'/e 
are significantly below the data. This is an important result as it shows that even if the 
value of B^ 2 ^ from lattice calculations would move up in the future, the SM would face 
difficulty in reproducing the data provided the large -N bound in (4) is respected. 

With these results at hand, we are in the position to summarise the present picture 
of the estimate of e'/e in the SM: 

• First, parametric uncertainties decreased by much since the analyses of e'/e around 
the year 2000. This includes the uncertainty in ImA^ which is presently about ± 7% 


and is irrelevant in the estimate in (61) but plays some role when e'/e is larger. 


Also the improvement on m s should be appreciated, entailing that the uncertainty 
on m s no longer is an issue. 

( 3 / 2 ) 

Second, the previously sizeable uncertainty due to B$ 7 7 has become sub-dominant, 
much smaller for example than the one due to isospin violation. This is thanks to 
impressive progress on the lattice 23 , which confirms large -N estimates employed 


in our previous papers, but with far smaller uncertainty. 
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• Third, the present analysis further increased the effectiveness of our framework, 
leading to a situation in which a single parameter B^ 2 ' 1 is playing the decisive role 
in the answer to the question whether s' je in the SM can be reconciled with the data 
or not. The new finding both by the lattice QCD and large-IV approach that B g 1 ^ 
is below unity narrowed significantly the range for s' je in the SM in our framework. 

This picture clearly indicates the emergence of a new anomaly in K physics. As this 
anomaly is strictly correlated in our framework with the value of B^ 2 \ this parameter 
must be a priority for future non-perturbative calculations for flavour physics. Fortu¬ 
nately, it is accessible by first-principle lattice-QCD calculations. Systematic improve¬ 
ment is hence possible. (See also comparison with lattice below.) Progress on isospin 
violation will also be important. 

But already now, the results presented here motivate further scrutiny of the SM pre¬ 
diction as well as searching for viable beyond-SM explanations. We will briefly discuss 
both directions, in Sections [5] and [6, respectively. 

Last but not least, the great reduction in parametric and hadronic uncertainties, made 
effective through our formalism, and good prospects on B^ 2 \ may make a more precise 
measurement of e'/e in the future worthwhile. 


3.2 Discussion of Bq 1 ^ 2 ^ dependence 

The domination of our error estimate by the uncertainty on B g 1 ^ 2 * leads us to investi¬ 
gate the dependence of e'fe on B^ 2 ^ in more detail. 

There is a hierarchy in the four contributions discussed in the previous section with ii) 
being most important followed by iv), i) and iii). For central values of input parameters 
we find 


- = 10 

£ 


-4 


ImAt 


1.4 • 10 


-4 


(1 - D eff ) ( - 4.1(8) + 24.7 B£ /2) ) + 1.2(1) - 10.4 B 


?(3/2) 


(64) 

with the four terms corresponding to the four contributions in question. The first number 
in brackets comprise the uncertainties of the sub-leading hadronic parameters q , p 3 , p 5 , 
p 70 and B^ 2 \ while the second number in brackets is due to the uncertainty in p 72 . This 
assignment of uncertainties will simplify the comparison with (68), even though it does 
not strictly follow our formalism. Furthermore, a remark on error correlations is in order. 
Due to implementing the constraints from CP-conserving data, correlations between the 
different contributions to e'/e are introduced. However, as the initial correlations of the 
hadronic matrix elements determined on the lattice are not available, we refrain from 
incorporating them into our analysis. 

It should be noted that the term representing Q e penguin operator involves the product 
a(l — Vt e s)B^ 2 \ Therefore, effectively isospin breaking corrections lower the value of 
B g 1 ^ by 0.866, implying in the case of B g 1 ^ = 0.57 an effective value of 0.49. 

In Figure [lj we show e'/e as a function of B g 1 ^ 2 '* for different values of B g 3 ^: 


b£ /2) = 0.7 (blue), Bf /2) = 0.8 (red), 
Bg 3//2) = 0.9 (green), B^ 2) = 1 (brown) . 


(65) 

( 66 ) 
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Figure 1: e'/e as a function of B^ 2 \ For further explanation see the text. 


The vertical band represents central value and error on from (3), the horizontal 

band the experimental world average on s' je. The black region on each line is excluded 
by the bound (4). We observe that the experimental value of e'/e can only be reproduced 
in the SM far outside the RBC-UKQCD range and then only for values B g 1 ^ > B g 3 ^ 
and B g 1 ^ > 1 in variance with the bound (4). 

/i/o'i ( 3 / 2 ) 

We finally observe that even if the bound Bq < B$ is violated, but the bound 
B^ 2 ' 1 < 1 is respected, the SM cannot quite reach the experimental data. Indeed, em¬ 
ploying this unlikely hypothesis, we find this time 


(eVe)sM = (H.l ± 3.2) x 10" 4 , ( B ( 6 1/2) = 1.0, Bf /2) = 0.76). 

4 Comparison with RBC-UKQCD lattice QCD 

4.1 Preliminaries 


(67) 


(64), 


The results for e!/e presented in 23,25 can be summarised by a formula analogous to 

ImAt 


- = 10 

£ 


-4 


1.4 • 10 


-4 


6.5(3.2) + 25.3 B {2,2) + 1.2(8) - 10.2 B, 


(3/2) 


( 68 ) 


In deriving this formula, we used the value of the matrix element (Qq) given in 25 
for /j, = 1.53 GeV: 


(QM)o = - 0.379(97) (83) GeV 3 (RBC-UKQCD) 


(69) 






















4 Comparison with RBC-UKQCD lattice QCD 


18 


with the first error being statistical and the second systematic. Using (41), we find (see 
also (24]) 


B, 


( 1 / 2 ) 


(/i = 1.53 GeV) = 0.57 ±0.19, 


(70) 


and consequently (3). The value of B g 1 ^ 2 * is significantly lower than its upper limit from 

0 |24] and the values for B'^ 22> used in many papers until now. 


U m 


the large -N approac 
This is the central reason why the lattice result is substantially below the data. 


Using (42) and comparing to the corresponding matrix element in 123 one extracts |21| 
Bf ,2 \3 GeV) = 0.75 ± 0.05, Bf /2 \m c ) = 0.76 ± 0.05, (71) 


which displays the very weak /j dependence mentioned above. 

Setting B^ 2) = 0.57 ± 0.19 and Bf^ 2) = 0.76 ± 0.05, we indeed obtain the result i 

©>• 


m 


Comparing formulae (|64) and (68), we observe the following differences: 

In 


25 


, a = 1 and O e ff = 0 have been employed. 


The main difference for fixed B^ 2 ' 1 and B^ i, " 2 ' > is found in the first term in (|68 


Not 


only is the error in this term much larger than in our formula but also is this term 
significantly larger than found by us. 

• Also striking is the sizeable error in the third term which is very small in our case. 

Let us then have a closer look at the contribution of the Q 4 operator in order to clarify 
the reason for this difference. 


4.2 Contribution of Q 4 and ReA 0 

Using the formulae of the previous section, we readily find 

2 y /2 ReA 0 


{Qi(m c )) 0 = . . 

(1 ± q)z _ GpV U( iV us 

For q = 0.05, using the experimental value of Redo; we obtain 

(Qi(m c )) 0 = 0.22(1) GeV 3 . 


On the other hand in 25 q ~ 0 and 

Re.4 0 = 4.62(0.95)(0.27) x 10“ 7 GeV, 


(72) 


(73) 


(74) 


the central value of which is roughly 40% larger than the experimental value in (47). From 

(75) 


(72) we now find 


(Q 4 K))o = 0.31(7) GeV 3 . 

This value agrees with the one given in | |25] : 

(Q 4 (1.53 GeV)) 0 = 0.271(93)(60) GeV 3 (RBC-UKQCD). 


(76) 
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But what is striking is the high precision obtained for this matrix element in our 
approach and still large uncertainty in the lattice result. It should also be noted that 
the contribution of the operator to e'/e is in the present lattice result comparable to 
the one of Q 8 and can be even larger than the latter one, which is not possible in our 
approach. 

4.3 Electroweak contribution 


On the other hand the electroweak penguin contri 
lattice value for ReA 2 agrees well with experiment 
0.76 ± 0.05, we find 


nution to e'/e is similar because the 
23]. Using the lattice result B g 3 ^ = 


(e'/e) E wp = - (6-7 ± 0.5) x 10 ” 4 , 


(77) 


which can also be obtained from the last two terms in (64). This result compares well 
(e'/e) E wp = - ( 6.6 ± 1.0) x 10" 4 , (RBC-UKQCD), (78) 


with 23 


although our error is substantially smaller. 


5 Can the large observed e'/e be made consistent 
with the SM? 

Given the significant discrepancy between our SM prediction and the experimental result, 
we first consider possible missing or underestimated contributions in the SM. 


5.1 Missing chromomagnetic contributions 


In our discussion (and much of the literature) the chromomagnetic penguin Q$ g (and also 
its electromagnetic counterpart have been tacitly dropped. It is straightforward to 
extend the formalism to include Q$ g , which being pure A/ = 1/2 impacts only on Inn4 0 - 
While ?/ 7 7 is small compared to the leading electroweak penguin coefficients, precluding 
any effect, the coefficient y% g is sizeable. The status of the hadronic matrix element 
( Qsg)o is rather uncertain. A calculation at leading non-vanishing order in the chiral 
quark model 39 gave 


(Q8g)o — — h 


1 11 


77U 


Fk 


16 7 r 2 2 rn s + 


m K rn, 


B : 


8 g 1 


(79) 


(recall h = yj 3/2 in our normalisation) with B$ g = 1, obtaining an upward shift of about 
0.3 x 10 ~ 4 on e'/e. Due to uncertainties from unknown higher orders and 1/A corrections, 
an ad-hoc range 1 < B Sg < 4 was advocated in 40 for setting bounds on new physics. 
For Cs g (m c ) ~ —0.185 and central values of our other input parameters, the resultant 
shift is in the range 


a£ 

£ 


= (0.2... 0.7) x 10 


-4 


(80) 


QSg 
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At the upper end of the range, while still being insufficient to explain the tension between 
theory and experiment, the contribution becomes competitive with some of the larger sub¬ 
leading uncertainties. Although a chromomagnetic contribution has never been seriously 
considered as a sizable SM contribution, the possibility cannot be fully excluded. A more 
definite conclusion would be desirable and will require the computation of ( Qs g )o hi the 
large-At approach or on the lattice. 


5.2 Missing low-energy contributions 


Two of the largest terms in our error budget concern low-energy physics: hadronic matrix 
elements in the isospin limit, as well as corrections to the isospin limit. In 41,42 it has 


been pointed out that for approaches that do not include final-state interactions, analytic- 
ity suggests extra positive contributions to the value of B^ 2 ' 1 and negative corrections to 

(o /o') 

the value of B% 7 , both of which would raise e' je. (See however 
the correction factors of 


41 


42 


43 


If we naively apply 

to typical large-A" values B^ 2] = 0.6 and B^ 2 ' = 0.8, 
an increase of e' je to 7.8 x 10~ 4 results, still well below the data. (Employing the lattice- 
inspired central values in our error estimate, B^ 1] = 0.57 and B^~ ] = 0.76, results in 
a very similar value e'/e = 8.5 x ICE 4 .) While a complete non-perturbative calculation 
should account for the full matrix elements including final-state interactions, the issue of 
final-state interactions may not yet be completely under control 3 and certainly deserves 
further study. 

Another type of long-distance corrections is isospin breaking, both due to electromag¬ 
netism and rn u ^ rrid■ This is parametrised by the two parameters f2 e fr and a. The latter 
only affects the overall normalisation and cannot bring the SM into agreement with data. 
Explaining the measured e' je due to the former would require a value of opposite sign 
and an order of magnitude larger than the value obtained in 29,30 . Nevertheless, given 


the profound implications of the e' je anomaly, this issue deserves further scrutiny, and 
also lattice-QCD studies should take these corrections into account. 


5.3 Missing higher-order corrections to the Wilson coefficients 


Higher-order corrections to the Wilson coefficients will also have an impact on the theory 
prediction of e' je. While it seems highly unlikely that they can bring the SM prediction 
into agreement with experiment, it is still instructive to discuss them in slightly more 
detail. In our analyses we fixed the renormalisation scale to /i = m c in the three-flavour 
theory. Hence the computation of the Wilson coefficients involves several steps, which 
start with matching at the weak scale and end with integrating out the charm quark at 
H = m c . The intermediate steps involve the renormalisation group evolution of Q\ - Qio 
and integrating out the bottom quark. 

The weak-scale matching corrections are known at NNLO for the electroweak pen- 
as well as the current-current and QCD penguins 


12 


44 , albeit in a different 


gum 

renormalisation scheme for the later two. The respective scheme transformation is given 
13 , where the relevant anomalous dimensions for the NNLO evolution of Q\ - Qq can 


m 


d For instance, the final-state phase shifts obtained in 25 are not in good agreement with the values 


extracted from experiment. We thank Chris Saclirajda for discussion. 
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also be found. For these operators the matching corrections at fj, = nib are also known 14 


yet all other matching corrections and anomalous dimension matrices are currently known 
only at NLO. 

In particular the unknown matching corrections at /i = m c could be sizeable 38 


since the strong coupling is growing rapidly in this region. For this reason we estimated 
higher-order corrections by varying the matching scale around /i = m c , and used the three- 
flavour renormalisation group running to determine the Wilson coefficients at // = m c . 
The resulting residual scale dependence is typically in the ball park of 10% for 2/3-7/10, 
but substantially smaller for z + and Z-. Using this procedure, only the uncertainties in 
2/6, and to a lesser extent y 8 , have a significant impact on the error budget of e'/e. 

The partially known NNLO corrections to y 8 are quite large 


12 


and decrease the SM 

prediction for e'/e. Accordingly, only the NNLO corrections to y 8 could arguably lead to 
a significant enhancement of e'/e, but our error estimate shows that a 10% increase in y$ 
results only in a 1.2 x 10 - ' 1 increase in the SM prediction. Bringing the SM prediction close 
to the experimental value would require a very large higher-order correction to y$ which 
would cast serious doubts on the convergence of the perturbative series in our approach. 
If this was indeed the case, we would have to perform our analysis in a four-flavour setup, 
i.e. above the charm scale, which would also require new calculations of matrix elements 
on the lattice. 


6 BSM physics in e'/e 

Not having been able to identify a plausible way to reconcile our prediction with the 
data (other than attributing it to a large statistical fluctuation somewhere), we turn to a 


discussion of physics Beyond the Standard Model (BSM) in e'/e. We first note that (19), 
reproduced here for convenience: 


£ 

£ 


UJ + 


y/2 \£k\ . 


Im/h 


(1 — fl e ff)- 


1 hn/U 


a ReA 2 


(81) 


remains intact in the presence of new physics, which can be classified by which of ReA 0j2 
and ImA 0;2 is affected. 


6.1 BSM physics in ReA 0 2 


Noting that the RBC-UKQCD prediction 25 of ReA 0 exceeds the experimental determi¬ 
nation, while the large-fV method exhibits a deficit 17 , we define the ratio 


H = 


(Re A 


OjSM 


(ReA, 


O/EXP 


(82) 


which takes the central value H = 1.4 and H = 0.7 in 25 and 17 , respectively. I11 


other words, we are considering a scenario where the experimental value of ReAo is a 
sum of the SM contribution and a BSM contribution. We cannot presently exclude that 
such a sub-leading part of ReA 0 comes from NP, a possibility investigated in 20 . As we 


have seen there is a strong correlation between ReA 0 and the matrix element of Q 4 and 
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Figure 2: e'/e as a function of Bg 1 ^, for three values of H defined in the text. 


consequently there is an effect on e'/e. We stress that the denominators in (19) are always 


the true (experimental) values including any BSM contributions. It is the numerator term 
ImAo that is affected through the correlation of hadronic matrix elements. 

Our formalism can easily be adapted to this case; one merely needs to multiply the 


V — A term given in (36) by a factor of H. In this fashion, the denominators in the ratios 
(36) are corrected for their BSM “contamination” and the theoretical SM expressions are 
again valid. Note that the ratio (38) is not modified. In Figure d, we plot e'/e as a function 
of B^ 2) for H = 0.7 (blue), H = 1.0 (black), and H = 1.4 (red). We see that taking 
the RBC-UKQCD central value for ReAo to be the true SM prediction, the agreement 
between theory and data for e'/e is worsened - and compensating for this requires even 
larger values of Bg 1 ^ 2 '* than in the SM. Conversely, taking the large -N central value at 
face value one observes a slight improvement (reduction) of the tension in e'/e by means 
of an upward shift. But in both cases, the effect is not huge, dwarfed by the uncertainty 
in B^ 2 \ and reconciling theory and experiment still requires Bg 1//2 ^ > 1. We conclude 
that CP-conserving data does not favour a scenario of BSM in Re/lo• although there is 
sizable room for it. A similar discussion could be given for NP in ReTh- 


6.2 BSM physics in ImA 0 ,2 

The result obtained in our paper that e'/e in the SM is significantly below the ex¬ 
perimental data has an impact on various NP models. This is in particular the case for 
models in which there is a strong correlation between e'/e and the branching ratios for 


rare decays K + —>• n + vv and K L —> 7 t 0 z/za Such a correlation has been stressed first in 45 
and investigated in many papers since then. See 46 and references to earlier literature 
therein. 
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In several models, like littlest Higgs model with T-parity (LHT) [47j, and generally 
Z-models with new FCNCs, only in left-handed currents [20^48], enhancement of the 
branching ratio for Kl —y is significantly constrained by e'/e because in these 

models such an enhancement is correlated with the suppression of e'/e with respect to 
the SM. This is also the case of K + — y 7 t + uu but as K + —y 7i + uu receives in addition to 
imaginary parts of the relevant amplitudes also the real parts, this correlation is much less 
pronounced. Therefore in such models in order to have large enhancements of Kl —> tt°uu 
and K + — » ■k+vv the SM prediction for s’ je must be above the data, which is certainly 
not favoured by our analysis. 

Therefore, in these models the agreement with the data for e'/e can generally be 
obtained only with strongly suppressed branching ratio for Kl —> n°uu. In the case of 
K + —> 7 t + uu this suppression is not required but significant departures from the SM are 
not allowed. The recent analysis within the LHT model 49 shows this explicitly. 

Now, in the models just described NP enters s' je only through ImA 2 and the presence 
of only left-handed FCNCs implies uniquely the strict correlation between e'/e and Kl —> 
tt°uu mentioned above. But as shown in 50 in the presence of both left-handed and right¬ 


handed FCNCs it is possible to arrange these couplings without significant fine-tuning so 
that the enhancement of e'/e required to fit data implies automatically the enhancement 
of Kl — >• 7r°z/P and to lesser extent of K + — » ti + vu. An explicit example of a model with 
tree-level Z exchanges contributing to e'je, Kl —>• tt°uu and K + —y n + i , u can be found 


m 


50 


In Z' models the situation can in principle be different even if e' je is only modified 
through I111A2 because flavour diagonal quark couplings to Z' could have proper signs 
so that e' je and Kl —> can be simultaneously enhanced in models with only left- 

handed flavour violating Z' couplings. As pointed out in 51 some 331 models have this 
property. 

Another route towards the enhancement of e'/e, less studied in the literature, are Z' 
tree level exchanges with flavour universal diagonal couplings to quarks. In this case 
IU1A2 is not modified and NP enters only ImAo through QCD penguin contributions. As 
demonstrated in 20,50 also in this model e'/e and K —» nvu can be simultaneously 
enhanced. Moreover, this can be achieved with only left-handed FCNCs. If the e'/e 
anomaly will be confirmed and future data on rare decays will exhibit such enhancements, 
models of this kind and the ones mentioned in previous paragraph will be favoured. 

Clearly there are other possibilities involving new operators, like supersymmetric mod- 

left-right models [5], but this is another 


54 


els [40j|52||53l, Randall-Sundrum models 
story which requires further study. 

7 Summary and outlook 


or 


Motivated by the recent results on K —> rnr amplitudes from the RBC-UKQCD collab¬ 
oration, we gave another look to the ratio e'/e within the SM. The main result of our 
analysis is the identification of a possible new anomaly in flavour physics, this time in K 
physics. This was possible because: 

• Improved results for the parameters B^ 2 ' 1 and B g 3 ^ became available through 
recent lattice-QCD studies by the RBC-UKQCD collaboration that are supported 
by the large -N approach which provides upper bounds on these parameters. 
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• We employed a formalism that is manifestly independent of the values of leading 
(V — A) Cg) (V — A) QCD penguin, and EW penguin hadronic matrix elements of 
the operators Q 4 , Q 9, and Q 10. In this manner a prediction for e'/e could be made 
that is more precise than presently possible by direct lattice-QCD simulations. 

I11 this context, we have presented a new analytic formula for e'/e in terms of B^ ,l2} 
and B g 3 ^ that is valid also in SM extensions with the same operator structure. This 
formula depends on the Wilson coefficients of the contributing operators which are model 
dependent while B g 1 ^ and B^ 2 \ related to long-distance dynamics, are independent of 
NP contributions. Thus our formula can be used for models such as the models with 
constrained MFV, 3-3-1 models and littlest Higgs models. We have also provided an 
update of a formula for e'/e in which NP enters directly through the shifts in basic one- 
loop functions. 

Our analysis emphasises the correlation between the amplitude ReH 0 and the contri¬ 
bution of the Q 4 operator to e'/e given in (72). As the central value of ReAo hi 25 is by 


40% above the data, this calculation overestimates the contribution of Q 4 to e'/e making 
it smaller. Assuming that ReAo is fully described by SM dynamics, we could improve the 


accuracy of its estimate by roughly an order of magnitude, as seen in (73) and (76). 
We have extracted from 


25 


the value of B g 1 ^ obtained by the RBC-UKQCD col¬ 
laboration to fold that it is significantly lower than unity. In fact this is the main reason 
for the low value of e'/e found in that paper. On the other hand, should the values of 
B g 7 7 and 7 7 eventually turn out to be close to the upper bound from the large -N 
approach 24 , significantly larger values of e'/e are found, although still roughly by a 


factor of two below the data. 

Our improved anatomy of e'/e clearly demonstrates that the SM has potential difficul¬ 
ties in describing the data for e'/e. However, there are several open questions that have to 
be answered before one can be fully confident that NP is at work here. Answering them 
would also allow us to give a better estimate of the room left for particular NP models. 

Our analysis shows that the next most important issues that have to be clarified are 
as follows: 

• The value of B g 1 ^ should be determined with an accuracy of at least 10%. Fig. ll 

( 3 / 2 ) 

demonstrates this need clearly, but also a higher precision on B g 7 7 would be ben¬ 
eficial. 

• The values of the Wilson coefficients Hi at the NNLO level. First steps in this 


direction have been taken in 12,13 


Improved calculations of isospin breaking effects, represented in our formula in (19) 
by the parameters a and O e ff- 

The role of electromagnetic corrections to the hadronic matrix elements, as empha¬ 
sised already in 10 . Without these corrections there remains some uncertainty due 


to the renormalisation scheme used for operators. 

Precise theoretical predictions of ReAo and Re%2 within the SM, which would tell us 
to which degree our assumption of neglecting NP contributions in these amplitudes 
is justified. 
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Finally, our understanding of the role of final-state interactions in e'/e, see |4| 
references therein, should be improved. 


and 


Onr present results could be affected to some extent by the future finding that some 
part of the amplitude ReAo does not come from the SM dynamics but NP. As Fig. 2 
shows, if ReAo in the SM is below the experimental value, as indicated by the large-iV 
approach 17 , the suppression of e'/e by the Q 4 operator is smaller, implying a larger 
value of e'je. On the other hand if ReA 0 in the SM is above the experimental value 
as presently seen in lattice data, the role of Q A will be enhanced and consequently ef/e 
smaller. But as seen in the error budget of Table [4, this effect is by far less important 
than the sensitivity to B^ 2 \ 

In view of the tendency of e'/e in the SM to be significantly below the data, it is 
exciting that in the coming years LHC might tell us what this physics could be. But also 
independent studies of e'/e in various extensions of the SM could select those extensions of 
the SM in which e'/e could be enhanced over its SM value. In fact first phenomenological 
implications of our results on new physics models have been presented in 50,55 . In any 


case it appears that e'/e could soon become again a leading light in flavour physics. 
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A Subleading contributions to ImAo,2 and related 
operator matrix elements. 

Including operators with small Wilson coefficients and colour-suppressed hadronic matrix 
element, the isospin ratios (36) and (38) receive the following corrections: 

( Inn4 0 


A 


A 


\ RgAo J y_A 

( IhiAq 
\R eA c 


= Imr 


P3 ( 2 ( 2/3 + Vi) ~ b (y 9 + y w )) 


0 / 6 


— — —7— J-Hl A/ 

V2 


2(1 + q)z. 

P5 2/5 (Qe)o + b ( 2 /s + pro 2 / 7 ) (Qs)o 
ReAn 


where we have defined 


Ps = 


($ 3)0 

(Q -)0 ’ 


P5 = 


(<55)o 
(<5e)o ’ 


P 70 = 


{Q7) c 
(Qs)c 


(83) 

(84) 

(85) 

(86) 


All three ratios are formally at least 1/iV-suppressed and multiplied by small Wilson 
coefficients. Note that in (39) we have already included the i/ 7 (Q 7)2 contribution into r/g ff ; 
eq. (39) then does not receive an additional correction. We also define B ^ 2 ' 1 through 


(Qs(aO)o — 2 h 


m 


K 


_m s (/i) + m d (n) 


E 1 d(1/2) 
^ 7 T £>8 


(87) 


This convention deviates from (5.12) in [101 but is motivated by the result 24 


B< I/2) = 1 + of - ) =1.1±0.1, 


( 88 ) 


strongly supported by the lattice result for (Q s(/i))o in |25 


To keep our phenomenological formulae simple and central values transparent, we set 
the central values of p$, p§, and P 70 to zero and allow generous error ranges that comprise 
both the intervals expected from large-IV counting and those computed in 25 . For the 


ratio P 72 defined below (40), which is also 1/iV-suppressed and plays a very minor role 
numerically, we take the central value from |25] and, conservatively, attribute a 100% 
error to it. We furthermore employ B^ 2 ^ = 1.0 ± 0.2, also derived from 25 . All input 
rages are summarised in Table (2. In this treatment, we tend to overestimate our error 
on e'/e, but as shown in the body of the paper, this has a very minor impact on our 
predictions. 


B Analytic Formula for e r /e 

The expression ([50) can be put into a formula that is more useful for numerical evaluations 
as it shows explicitly the dependence on rrit and m s . The most recent version of it has 


been presented in 20 , but we update it here due to the change of some input parameters 
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entering the formulae for hadronic matrix elements and a different treatment of isospin 
breaking corrections. We then have 


= IrriA t ■ F e '(x t ), 


(89) 


SM 


where 

F £ \x t ) = P 0 + Px X 0 (x t ) + Py Y 0 (x t ) + P z Z 0 (x t ) + Pe E 0 (x t ) ■ (90) 


The first term is dominated by QCD-penguin contributions, the next three terms by 
electroweak penguin contributions. The last term expresses the m t dependence from con¬ 
tribution of QCD penguin operators and is totally negligible. The x t dependent functions 
are given as follows 


X 0 (x t ) 


Y 0 (x t ) 


Xt 

8 

Xt, 

8 


x t + 2 

x t - 1 

x t — 4 

-7 + 

x t - 1 
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3 x t 
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In x t 


(91) 

(92) 


Z 0 (x t )= - - In x t + 


1. 18xf - 163x 3 + 259x 2 - 108x t 


144 (x t - l) 3 


+ 


32xf - 38x 3 - 15x? + 18x t , 

+ 72 (x t - l) 4 n ^’ 

2 x 2 t (15-16x t + 4:x 2 t ) x t (18 - llx t - x 2 ) 

W = - 3 lnXt + -6(TV7]4- ln x ‘ + 12(1-x,)3 

where x t = m 2 /M 

The coefficients Pj are given by 

P^rf + rfPe + rfPg, 

where we have defined 


— d(1/2) 


P-6 = Bq 


114.54 MeV 
_m s (m c ) + m d (m c )_ 


— d(3/2) 


Ps = B 


114.54 MeV 
_m s (m c ) + m d (m c )_ 


(93) 

(94) 

(95) 

(96) 

(97) 


The coefficients and comprise information on the Wilson-coefficient func¬ 

tions of the AS = 1 weak effective Hamiltonian at NLO and incorporate the values of the 
matrix elements of those operators that we could extract by imposing the experimental 
values of Re/lo and Re/l 2 • Their numerical values are given in the NDR-MS renormalisa¬ 
tion scheme for fi = m c and three values of a s (M z ) in Table [5| 
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a s (M z ) = 0.1179 

a s {M z ) = 0.1185 

a s {M z ) = 0.1191 
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-1.788 

0.665 
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Table 5: The coefficients r[°\ r- 6) and r- 8 ' 1 
three values of a s (M z ). 


of formula (95) in the ADR-MS scheme for 
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